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Abstract 

While the standard (introductory physics) way of computing the equvalent resistance of non¬ 
trivial electrical ciruits is based on Kirchhoff’s rules, there is a mathematically and conceptually 
simpler approach, called the method of nodal potentials, whose basic variables are the values of 
electric potential at the circuit’s nodes. In this paper, we review the method of nodal potentials and 
illustrate it using the Wheatstone bridge as an example. At the end, we derive—in a closed form— 
the equivalent resistance of a generic circuit, which we apply to a few sample circuits. The final 
result unveils a curious interplay between electrical circuits, matrix algebra, and graph theory and 
its applications to computer science. The paper is written at a level accessible by undergraduate 
students who are familiar with matrix arithmetic. For the more inquisitive reader, additional proofs 
and technical details are provided in the appendix. 
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I. INTRODUCTION 


The first thorough mathematical description of electrical circuits goes back to Gustav 
Kirchhoff.^ Ever since, this topic has given rise to a great number of pedagogical arti¬ 
cles. On one hand, circuits are a source of elegant problems and solutiond^ or interesting 
experiments.®^ On the other hand, electrical circuits can be used to visually and intu¬ 
itively illustrate some complicated (or often misunderstood by student^^ electrodynamical 
concepts.^ Some authors apply powerful methods known from elsewhere to compute the 
equivalent resistance of nontrivial electrical circuits.^®®^ In this paper, we primarily focus 
on computing the equivalent resistance of a generic circuit. Remarkably, unlike many of the 
works cited above, the method we employ here doesn’t require the individual resistances to 
have specihc values, e.g., to be all equal, or to form a symmetric or lattice-like structure. At 
the same time, this method is simple enough to be accessible to students familiar with the 
basics of matrix arithmetic. 

From the mathematical perspective, a resistive electric circuit can be understood as 
a graph whose nodes are connected via links (edges) with assigned numerical values of 
resistance (R). In addition, two nodes are assumed to be “connected to the battery,” which 
hxes the potential difference (voltage, £) between the two nodes. The battery gives rise to 
currents flowing along the edges of the circuit. These edge currents (/) are related to the 
potential differences across the corresponding nodes (AU) via Ohm’s law: 


I = 


AU 


( 1 ) 


The total current flowing through all the edges coming out of a battery node is called the 
output current (/out)- The total current that flows into the other battery node is equal to the 
output current. One of the important characteristics of an electrical circuit is its equivalent 
resistance (/?eq)- If is defined as the resistance of a single resistor (edge) that, if it were to 
replace the whole circuit, would result in the same amount of the output current. In other 
words, 

Rcq = J • (2) 

-^out 

Although the output current is explicitly present in the definition of equivalent resistance, 
it is clear (e.g., on dimensional grounds) that Req depends only on the given edge resistances. 
In some simple cases, it can be computed without finding the output current. For instance. 
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if a node is only shared by two edges, such edges are said to be connected in series, and 
the individual resistances are merely added. If two resistors are on edges connected across 
the same pair of nodes such a connection is called in parallel, and the equivalent resistance 
is computed as the reciprocal of the sum of reciprocals. Both of these situations can be 
generalized straightforwardly for more than two resistors: 


i?eq = Ri + R 2 + ■ ■ ■) for connection in series; 

i?eq = (-Rr^ + ^ 2 ^ + • • •) ^, for connection in parallel. 


( 3 ) 

( 4 ) 


Some larger circuits allow for reduction by identifying sub-circuits whose elements are con¬ 
nected either in series or in parallel. Replacing such sub-circuits by single equivalent resis¬ 
tances turn by turn may result in a trivial circuit, thus giving an algorithm for computing 
equivalent resistance. 

At the same time, it is clear that not every circuit can be simplihed this way. For example, 
a circuit with no parallel edges and whose nodes are at least three-valent (such nodes are 
referred to as junctions) cannot be reduced in the above sense. The simplest non-simplihable 
circuit is depicted in Fig. It is easy to see that there are no elements connected in series 



FIG. 1: The simplest non-simplihable circuit (Wheatstone bridge). There are no 

connections in series or in parallel. 


or parallel. In order to determine the equivalent resistance of such a circuit, one typically 
introduces the unknown edge currents and writes down Kirchhoff’s rules. After solving those 
equations for the unknown currents, one can compute the output current and obtain the 
equivalent resistance from Eq. ([^. 

There is also an alternative approach to Ending unknown quantities associated with 
electrical circuits, called the method of nodal potentials (see, for example. Refs. |13| and | 
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there are also variational alternatives to Kirchhoff’s loop rule presented in Ref. [T5j). In 
this method the basic variables are the potentials at the circuit’s nodes rather than the 
edge currents. The advantage of this method is that it usually involves signihcantly fewer 
variables than in Kirchhoff’s framework. In addition, Kirchhoff’s rules come in two types 
(junction and loop) which have distinct algebraic structure. On the contrary, the equations 
for the unknown nodal potentials are all of one type (following from the junction rule). This 
presents a major simplihcation for the analytical description of electrical circuits, as well as 
for hnding various quantities of interest. 


This paper is organized as follows. We start by reviewing the method of nodal potentials 
for electrical circuits, highlighting some immediate implications of the formalism and working 
out a sample circuit (the Wheatstone bridge). For generic values of the edge resistances, the 
latter is somewhat nontrivial when using Kirchhoff’s rules,^but is quite straightforward us¬ 


ing the method of nodal potentials (see Eq. (17) for the equivalent conductance). In Section 


IIIB, we then proceed to derive a closed formula for the equivalent resistance/conductance 


of an arbitrary resistive circuit (see Eq. (25)). Related results exist in the graph theoretical 
literaturej ^^*^^ ! expressing resistance distance of a graph with l-fl resistors on each edge. Ref. 
[HI utilizes the Moore-Penrose (pseudo) matrix inverse and does not lead to easily derived 
conclusions. In this paper, the derivations, as well as the hnal result, are quite intuitive 
and accessible to undergraduate students familiar with basic matrix arithmetic. In the last 
section, we discuss generalizations of the formula, as well as its possible relation to matrix 
algebra and graph theory and its applications to computer science. 


II. REVIEW OF THE METHOD OF NODAL POTENTIALS 


In this section we review a convenient description of resistive DC-circuits containing one 
battery and investigate its implications. We shall see how the method of nodal potentials 
can simplify the analysis of an electric circuit compared to the standard Kirchhoff’s rules 
technique. The same approach can be applied to DC-circuits with several batteries, as well 
as to generic AC-circuits. Such generalizations are discussed in the last section. 
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A. Formalism and notation 


Consider a circuit containing n nodes, such that nodes 1 and n are connected to the 
positive and negative terminals of the battery respectively (Fig. [^. 
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FIG. 2: A generic circuit with n nodes. The positive (higher potential) battery terminal is 
connected to node 1 and delivers the output current to the circuit (/out)- The current flows 
back into the other terminal of the battery via node n. 


Every link connecting two nodes i and j is assumed to have known resistance Rij = Rji. 
In fact from now on, it will be more convenient to use the conductance rather than resistance, 
dehned as 

1 

( 5 ) 

rtij 

thus giving rise to the conductance matrix aij. It is easy to see that the rules for computing 
equivalent conductance are reversed compared to those for equivalent resistance: 


o'eq = cTyi + (Ts + ..., for conuectiou in parallel 

(Teq = + ...) ^, for connection in series 

Without loss of generality, we can make the following assumptions: 


( 6 ) 

( 7 ) 


• Vi = and Vn = 0. Since electric potential is dehned up to a constant, we hx one of 
them to zero. The potential at the positive terminal, S, can be used as a unit. 

• Every node is connected to every other node. If, in reality, some nodes i, j etc. do 
not share a link, we simply put CTjj = 0. 

The latter assumption will allow us to not worry about the circuit’s topology and concentrate 
on purely algebraic description. In fact, if we only focus on non-simplihable circuits, every 
node, except possibly the 1®* and ones, will have at least three edges. Indeed, a node with 
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only two edges would imply a connection of the edges in series which could be replaced by 
a single edge. Note that each column/row of the conductance matrix for a non-simplihable 
circuit must have at least three non-zero entries. 

We now dehne the edge current lij as the current flowing from node i to node j. Its 
expression can be readily written in terms of the nodal potentials as 

hj ■= (^ij {Vi - Vj) = -Iji, ( 8 ) 

which is manifestly anti-symmetric. In particular, the “diagonal” current la = 0. Since 
electric current flows from a higher potential to a lower one, this dehnition sets an outgoing 
current to be positive, whereas an incoming current would be negative. 

We conclude this subsection with the following remark. A minimally connected, non- 
simplihable circuit has at least ~ 3n/2 edges, hence the same number of unknown currents. 
At the same time, the number of unknown nodal potentials is (n — 2), which can be sub¬ 
stantially less than the number of currents. Thus the method of nodal potentials deals with 
fewer variables at the on-set. 


B. Implications 

Assuming that the variables describe a real circuit which has specihc unique values of the 
nodal potentials, Kirchhoff’s loop rules will be automatically satished by construction: any 
closed loop will come back to the same value of potential, hence making the overall potential 
difference zero. 

For a generic circuit, Kirchhoff’s junction rules take the form 

n n 

lij = '^{Vi- Vj) (Tij = 0, for i 7^ 1 or n. (9) 

i=i i=i 

Here, each sum constitutes the total nodal current for i = 2, ...,n — 1. Nodes 1 and n are 
not included, as the battery provides non-zero nodal currents at those nodes. In fact, Eqs. 
(|^ imply the following 

Lemma 1 The total current flowing into the node is equal to the total current flowing 
out of the 1®* node. 
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Proof: Since the edge currents are anti-symmetric, lij = —Iji, the sum over both indices 

n 

'Yf hj = 0. Splitting the summation over i into i = 1, i = n and the rest yields 
*,i=i 

n n n—1 ( \ 

-^1 j+x] ^n,j +f I=o- (lo) 

j=l j=l i=2 \j=l J 

By virtue of Eq.0 , each summand in the last term vanishes. The hrst and second terms 
are the output and (negative) input currents respectively, which proves the lemma. □ 

While the statement of the lemma was not mathematically obvious from the construction, 
it makes clear physical sense: since there is no accumulation of charge in the circuit, the 
total current coming out of one terminal of the battery has to equal the current flowing back 
into its other terminal. 

There is another statement which is obvious from the physical point of view, but non¬ 
trivial mathematically® 

Lemma 2 The values of the nodal potentials in a resistive circuit connected to a single 
battery should lie (strictly) between the lower battery voltage and the higher one, i.e., between 
0 and £. In other words, the battery sets the lowest and the highest possible potential in the 
circuit. 

Proof: Suppose that the maximum potential is attained at node m ^ 1. Then all the 
neighboring nodes would have a lower potential resulting in currents from node m 

^mj ~ CTmj (hm ^ 0 (H) 

to be outgoing. The latter would violate the junction rule (|^. Thus the maximum potential 
can’t be attained at any node other than the 1*^* one. Similarly, the minimum potential is 
attained at node n. Hence all the nodal potential values lie between Vi and Vn, which proves 
the lemma. □ 

In the next session we illustrate how the method of nodal potentials helps to determine the 
equivalent conductance of the Wheatstone bridge circuit with arbitrary edge conductances. 
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C. An example: Wheatstone bridge circuit 


We now assume that the values of edge resistance in Figjljare given and hnd the equivalent 
resistance (conductance) of the circuit. Notice that if one was to solve this problem using 
Kirchhoff’s rules, one would have to deal with six unknown currents: one in each resistor 
plus the output current. At the same time, there are only two unknown nodal potentials: 
V 2 and V 3 , since Vx = and V 4 = 0. For convenience, we also label the edge conductances 
similarly to the original resistances: 

_ 1 _ 1 _ 1 _ 1 _ 1 
<^12 = <^24 = cr 2 — —, CT 23 = CT 3 — —, cri 3 = CT 4 — —, CT 34 = CT 5 — —. 

ril ri 2 4 x 3 r£4 rts 

( 12 ) 

It is easy to see that there are exactly two junction equations for this circuit, one for 
node 2 and one for node 3: 


.^21 + f^24 + f^23 — O'! (14 — Vi) + (72 (V 2 — V 4 ) + (T 3 (V 2 — V 3 ) — 0 (13) 

Izi + Im + -^32 = <74 (V 3 — Vi) + (Ts (V 3 — V 4 ) + <73 (V 3 — V 2 ) = 0 (14) 

Setting Vi = ^ and 1/4 = 0 and solving the equations on the righthand side for the 

unknown potentials V 2 and V 3 , we obtain: 


_ f,0'lC3 + <73(74 _ j,<74C2 + (7i(73 

V2 — O - v^ — c - 

C2C3 - C2C3 - 

where C2 = cti + (72 + (73 and C3 = 0-3 + 0-4 + 


(15) 


It is important to understand that for a connected circuit (0 < aij < 00 ) the denominator 


m 


(15) can never be zero. To make it more transparent, we rewrite the denominator as 


T’ = cr3((7i + £72 + £74 + a^) + ((7i + C72)((74 + ( 75 ). 


(16) 


Since each term in T> is non-negative, it can only be zero if each term is zero. Irrespective 
of whether (73 = 0, the latter implies that ( 7 i = (72 = (74 = (75 = 0. This corresponds to 
a disconnected circuit, such that nodes 2 and 3 are completely isolated, which makes their 
potentials undetermined.^^ 

In order to hnd the equivalent conductance, it is easiest to consider the current howing 
into node 4, I 24 , + = <^2 (V2 — V4) -|- (V3 — V4). Setting V4 = 0 and using the nodal 





potentials in (15), we obtain 


<^eq — 


(T2V2 + 0-51/3 


0-lCr2C3 + a2<J3<Ji + 0-iCr3(T5 + 0 - 4(7502 


( 17 ) 


£ 0-3(0-! + 0-2 + 0-4 + 0-5) + (0-1 + 0-2) (0-4 + 0-5) ■ 

Looking at the answer, we see that it is a ratio of two polynomials of degree (n — 1) and 
(n — 2) respectively. This clearly guarantees the correct units of conductance. Moreover, 
each polynomial is a sum of non-negative terms. As explained in Ref.^ , this means that 
(15) and ([T^ are neither zero nor inhnity for any connected circuit. 

Finally, the Wheatstone bridge has a well known feature that for a special arrangement 
of edge resistances (conductances), there is no current through the middle wire (labeled with 
( 73 ). We can arrive this condition by setting equal potentials at the ends of the wire 14 = Va- 
This yields ( 74(75 = ( 72(74 or the standard 


R\ 

R2 R5 


( 18 ) 


In the next section we shall generalize the expression for the nodal potentials (15) and derive 
the equivalent conductance for an arbitrary circuit. 


III. METHOD OF NODAL POTENTIALS IN A GENERIC CIRCUIT 

We now revisit the generic circuit displayed in Fig. |^with the same assumptions as in 
Sec. IIA the battery terminals read Vi = £ and 14 = 0 and all the edge conductances aij 


are given. 


A. Expressions for the nodal potentials 


We can write the junction equations, analogous to (13) and (14), for all the nodal poten¬ 


tials, including Vi and 14- Collecting similar terms, we cast those equations in the following 
matrix form 


' Cl 

— (^12 

—O'lS 


— O'l^n 


' 14 ^ 


flout! 

— (^21 

C2 

—(^23 

■ —0'2,n-l 

— 0'2,n 


V 2 


0 

— C -31 

— <^32 

C 3 


~Cr3,n 


V 3 

= 

0 


— (^n-1,2 

1 

1 

CO 

Cjj_i 

^n—l,n 


V„_1 


0 

V ^n,l 

—0'n,2 

(^n,3 • • 

^n,n—l 

Cn j 


V Hi ) 


Vlin/ 


( 19 ) 
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where the unknown quantities are highlighted in boldface, and the diagonal elements Cj = 

n 

<^ij- We denote this matrix S,^and will also need S', its upper-left sub-matrix (n — 1) x 
i=i 

(n - 1) 


' Cl 

— Cu 

—Cls 

■ ■ 1 


/ 


IT \ 

• — <72,n-l 






C2 

— <723 

—C21 

C2 

—C23 

•• —C2,n-1 





—C31 

— C32 

C3 

■ ■ —C3,n-1 

and E" = 

— <732 

C 3 

• — <73,n-l 







— <7n-l,3 • • 

Cn—1 J 


S' = 


( 20 ) 

the lower-right sub-matrix of S' of size (n — 2) x (n — 2). Note that the sum of all the 
elements of S is zero, which implies that the n equations in ( |I^ are not independent. This 
is clearly consistent with the fact that Jin = —lout- Hence from now on we skip the 
equation. 

In order to solve for the unknown nodal potentials V 2 , V 3 , ..., W-i, we rearrange the 


equations in (19) as follows. The hrst element of each row of S multiplies Vi = £. We carry 


this term to the righhand side of each equation in (19). In the hrst row, we also carry the 


unknown output current lout to the lefthand side. Now these equations take the form 


^—1 — <7i2 — <7i3 ... — 


( lout \ 


( -ci\ 

0 C2 —<723 • • • — <72,n-l 


V2 


<721 

0 —<732 C3 ... —as^n-i 


V3 

= 8 

<731 

\ 0 1,2 1,3 • • • ^n—1 j 


\yn-ij 


^<7n—1,1 / 


( 21 ) 


Applying Cramer’s rule to rows 2 through (n — 1), we obtain the following expressions for 
the nodal potentials 




de^ 

det E"' 


( 22 ) 


Here S'^' is the matrix obtained from E" by substituting its column by ((T 21 , CT 31 ,..., an-i,iY■ 
On physical grounds, for any connected circuit the determinant in the denominator should 
be non-zero. This, however, is not so obvious from the mathematical point of view. See 
Appendix [A| for more detail. 


10 















B. Derivation of the equivalent conductance of a generic circuit 


In principle, one can find the equivalent conductance of the circuit in Figj^ using the 


same method as in Sec. IIC Specifically, we could use the nodal potentials found in (22) to 
compute the output current, using ([^ for n = 1, and substitute in 

lout 




(23) 


There is, however, a more economical way to arrive at the expression for the equivalent 
conductance. Interestingly, it can be obtained in a closed form. Moreover the expression of 
(Jeq in terms of the individual conductances (which are assumed given) is universal and does 


not require prior finding nodal potentials. Applying Cramer’s rule to the first row of (21), 
we can find lout 


.font — S 


-Cl 

— Cl2 

— Cl 3 

^1,71—1 

0‘21 

C2 

— C23 

• —0'2,n-l 

Csi 

— C32 

C3 


Cn—1,1 

1 

3 

1 

1 

1 

CO 

^n—1 


-1 

—C 12 

— <713 

<7l,n—1 

0 

C2 

— <723 

• —<72,n-l 

0 

— C 32 

C 3 

• —<73,n-l 

0 

— Cn-1,2 

— <7n-l,3 • ■ 

Cn—1 


= S 


- det S' 

— det S" 


(24) 


Therefore, the equivalent conductance reads 


det S' 


detE"' 


(25) 


As before, the answer is a ratio of two polynomials of degree (n — 1) and (n — 2), which 


clearly has correct units. Each determinant is of the form similar to that of Sec. IIC They 
both are non-zero (positive) for any connected circuit. See Appendix for more detail. 


Eq.(25) constitutes the main result of this paper. In the following section we consider a few 


specific electrical circuits and use (25) to compute their equivalent conductance. We discuss 


general properties of (25) in the last section 
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IV. EXAMPLES 


Even though the application of (25) is straightforward for an arbitrary electrical circuit, it 


is still interesting to consider a few examples and make sure that the formula work correctly. 


In this section, we analyze relatively simple circuits. More general properties of (25) are 
discussed in the Appendices. 


A. Wheatstone bridge 


Now, as we have the closed formula (25), we can compute the equivalent conductance of 


the Wheatstone bridge directly. Using the conductance labels of Sec. |II C[ the corresponding 
matrix is 

/ Cl —(Ti —(74 0 \ 

— (Ji C2 —CTs —(72 

— (74 —(73 C 3 —(75 

\ 0 -£72 -CTs C4 y 


E = 


(26) 


where ci = ai + ( 74 , C 2 = ( 7 i + (72 + ( 73 , C 3 = (73 + (74 + ( 75 , and C 4 = (72 + ( 75 . Applying ( [^ , 
we obtain 


Cl —(7i —(74 

-(7i C2 -(73 
— (74 —(73 C 3 


^eq 

C2 —C73 
— O'S C 3 

which agrees with ([T^. 


0'lO'2C3 + (72(73(74 + (7i(73(75 + (74(75C2 

cr3(cri + cr2 + CT4 + CTs) + (cTi + cr2) (cT 4 + CTs) ’ 


(27) 


B. A three-node circuit 


In Appendix B 1, we prove that formula (25) yields the right result when one combines 


two resistors in series and rewrites the S-matrix accordingly. For the three-node circuit in 


Fig. however, we shall obtain the equivalent conductance directly. Before we do so, note 
that - based on Eqs. 0 and ^ - the equivalent conductance is 


O'eq — O'lS 


0'12 + <^23 


(28) 
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FIG. 3: Simple circuits. The first and last nodes are connected to the battery. 


The E-matrix for this circuit is 


E = 


Cl 

— <712 

— <7i3 

—C'21 

C2 

— <723 

— C'si 

— <732 

C 3 


Setting the symmetric entries equal and using (25), we get 


<^eq 


Cl —0'12 
-(^21 C2 


C 2 


0‘120'23 + C'23CTi 3 -|- Cri2CTi3 
<^12 + <^23 


which is equivalent to (28). 


(29) 


(30) 


C. Edges without current 


In the four-node circuit of Fig. the current between nodes 2 and 3 is zero. Thus this 
circuit is electrically equivalent to the previous one. In particular, the equivalent conductance 
should not depend on (J 23 . It is insightful, however, to see explicitly how (J 23 drops out of 
the hnal expression. The E-matrix reads 


Cl 

— <7i2 

0 

— <7i4 

— <7i2 

C2 

—<723 

— <724 

0 

— <723 

C3 

0 

— <7i4 

— <724 

0 

C4. 


(31) 


Here we have already set the symmetric entries equal. Note that C 3 = a 23 - When computing 
the relevant determinants, we can replace the second row with the sum of rows 2 and 3, and 
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then do the same with columns 2 and 3. After these manipulations, the numerator of (25) 
reads 


det E' = 


Cl 

—C 12 

0 








Cl 

— <7i2 

—C 12 

C 12 + <^24 

0 

— <723 







— <7i2 

<7i2 + <724 

0 

0 

<^23 





(32) 


and (723 can be factored out. Repeating the same steps for det E", we can see that the 


denominator of (25) is also proportional to (T 23 . Therefore the equivalent conductance does 
not depend on this edge conductance. Clearly, the same argument works for any circuit 
containing an edge that only has one connected node. 


D. The classic resistor cube problem 

The edges of the cube in Fig. [^are made of l-f2 resistors {a = 1 We need to hnd 

the equivalent resistance of this circuit if the battery terminals are connected across the 
main diagonal (nodes 1 and 8 ). 


6 8 



FIG. 4: The classic resistor cube problem. Nodes 1 and 8 are connected to the battery. 
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The S-matrix for this circuit is 


3 

-1 

-1 

-1 

0 

0 

0 

0 

-1 

3 

0 

0 

-1 

0 

-1 

0 

-1 

0 

3 

0 

0 

-1 

-1 

0 

-1 

0 

0 

3 

-1 

-1 

0 

0 

0 

-1 

0 

-1 

3 

0 

0 


0 

0 

-1 

-1 

0 

3 

0 


0 

-1 

-1 

0 

0 

0 

3 


0 

0 

0 

0 

-1 

-1 

-1 

3 


( 33 ) 


Computing the relevant determinants, we obtain the equivalent conductance 


_ detS' 

d^ 


384 

3 ^ 


6 

5 




( 34 ) 


or Req = 5/6 n, which is the well-known result. Note that while symmetry is crucial for the 


classical solution of this problem, Eq. ( 25 ) would still apply for arbitrary edge conductances 
of the cube. 


E. Complete graph 


Consider a circuit with N nodes, such that every node is connected to every other node by 
an edge of hxed resistance, say 1 12 {aij = 112 “^). One method of calculating the equivalent 
resistance of this circuit is presented inland the answer is i?eq = 2 /iV ( 12 ). Before applying 


( 25 ), we note that there is a simple physical solution. Assuming that the battery is connected 


to nodes 1 and N, all the other {N — 2 ) nodes are equivalent, hence their nodal potentials 
are equal: V2 = C3 = ... = Vn-i- We can thus merge the latter nodes into one (M) and 
obtain a very simple circuit: {N — 2 ) parallel edges from 1 to M, [N — 2 ) parallel edges 
from M to N and one direct edge from 1 to N. It then straightforward to obtain the above 
expression for R^q. 

Let us now consider the relevant determinants of the E-matrix. 
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A- 1 

-1 

-1 . 

. -1 


-1 

A- 1 

-1 . 

. -1 

det S' = 

-1 

-1 

A -1 .. 

. -1 


-1 

-1 

-1 . 

. A -1 


( 36 ) 


with (A^ — 1 ) elements in each row and column. Expanding this determinant by the elements 
of first row, we obtain one term proportional to det S" and {N — 2 ) other terms, which can 
be shown to be identical (denote A). In other words. 


det S' = (A - 1 ) det S" + (A - 2 )A, 


(36) 


where 



1 

1 

. -1 


1 

1 

. -1 

det S" = 

... 1 

... 1 

. -1 

, A = 

-1 A -1 .. 

. -1 


1 

1 

. A -1 


1 

1 

. A -1 


( 37 ) 


with (A — 2) elements in each row and column. At the same time, the determinant in 


Eq.( 35 ) can be computed by subtracting the second row from the first one 



A 

-A 

0 .. 

0 


-1 

A -1 

-1 . 

. -1 

det S' = 

-1 

-1 

1 

. -1 


-1 

-1 

-1 . 

. A -1 


= A(det S" +A), 


(38) 


where we expanded the determinant by the elements of the hrst row. Setting the righthand 


sides of Eqs. ( 36 ) and ( 38 ) equal, we can express A as 

A = -IdetE". 


( 39 ) 


Substituting this result back into ( 38 ) yields det S' = ydetS". Finally, applying ( 25 ), we 
obtain 

det S' A ,„_i, 

^ d^ = y 

which agrees with the equivalent resistance i?eq = 2 /A (hi). 


16 
















V. DISCUSSION 


We have derived a closed formula for the equivalent conductance of an arbitrary circuit. 
All one needs to know is the edge conductances cxjj which give rise to the S-matrix (weighted 


Laplacian matrix) dehned in (19). The equivalent conductance can then be computed as a 
ratio of two sub-determinants of S via Eq. (25). Interestingly, both det S' and det S" are 
/mear functions of each individual edge resistance (see Appendix |A2[). 


One important feature of (25) is that it respects the permutation symmetry. Indeed, 
relabeling any nodes other than the two connected to the battery terminals (1 and n) must 
not affect the equivalent conductance of the circuit. For example, switching labels i and j 
would merely result in a minus sign in front of det S' and det S" thus keeping the answer 
intact. 

On physical grounds, the equivalent conductance of an arbitrary circuit (with cxjj < 
cxd) can be zero (for a disconnected circuit), but never inhnite. One the other hand, the 


denominator in Eq. (25) can vanish if sufficiently many atj equal zero. This implies that if 
det S" = 0, the other determinant, det S', must vanish as well. Physically, this is the case 
when the circuit contains isolated clusters of nodes. Removing such clusters (and reducing 
the size of the S-matrix) renders the determinants non-zero. 

In this paper we focused only on resistive circuits with a single battery. However, the 
same analysis can be applied to circuits with multiple batteries. This can be done by 
simply incorporating the additional EMF’s into the nodal potential differences. Clearly, 
ideal batteries would not affect the equivalent resistance between hxed nodes. 

Generalization to capacitive circuits is also straightforward, as equivalent capacitance 
obeys the same rules e and ([^ as conductance. Thus the hnal formula ( [2^ can be 
understood in terms of capacitance as well. 


In addition, Eq.(25) will work for the equivalent impedance (admittance) of an AC-circuit. 
The only difference would be that admittance is a complex number and either det E' or det E" 
can be zero even for a connected AC-circuit. In fact, setting the determinants to zero one 
can determine the resonance frequencies of the circuit. 


Interestingly, Eq. (25) unveils a curious interplay between electrical circuits, matrix 
algebra, graph theory and its applications to computer science. Specihcally, there is a 
straightforward correspondence between electrical circuits and random walks on graphs,^ 


17 











including the concept of escape probability^ which is a direct analog of equivalent resistance. 


In addition, Eq. (25) can help to investigate the connectivity of generic graphs, which is 
done in, e.g., lousing spectral analysis. These interdisciplinary connections are particularly 
useful, as there is much physical intuition about electrical circuits that could give rise to 
some less obvious mathematical statements. 

To conclude we note the following. While the method of nodal potentials should probably 
have place in an E&M curriculum,!^ I would rather not recommend handing formula (25) out 
to students; that is unless the students derive it themselves! Giving the answer away could 
ruin the excitement of puzzle-solving that students usually associate with circuit analysis. 
On the other hand, this closed expression can help physics instructors when creating multiple 
versions of circuits problems. For they would easily be able to randomize the edge resistances 
for many students without any worries about checking their students’ answers. 
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Appendix A: Properties of the determinants 


1. Positivity of detS' and det S'' 


In general, determinants contain both positive and negative monomials. However, in 
Sec. |II G we saw that, after some cancellations, the determinants in the numerator and 
denominator of dl?] ) had only positive terms. In this appendix, we shall prove that the 
determinant of S' is always of this form, that is 


Lemma 3 Each term in det S' enters with a plus. 
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Proof: We shall proceed with the proof by induction in n. For example, for n = 3, 


S' = 


Cl 

— C'21 



(Al) 


where Ci = a 12 + cxis, C2 = a2i + <723 and cri2 = <721. Thus the determinant 


det S' — C 1 C 2 — (7^2 ~ <^12(723 + cri2(7i3 + (7i3(723. 


(A2) 


In addition, in Sec. 11 C we saw the statement of the Lemma to be true for n = A. 


Assume now that the statement holds for n = k. In order to prove that it also is true for 
n = k + 1 (S' is then fc x fc), it is sufficient to demonstrate that the coefficient in front of each 
aij, as it enters det S', is positive. Without loss of generality, we can focus on ai2 = o'2i- In 
the determinant 


det S' = 


Cl 

—cri2 

— (^13 ■■ 


— C‘21 

C2 

— ^23 ■ ■ 

■ —0'2,k 

— 7‘31 

— 0'32 

C3 . . 

■ —<^3,k 

— <7A:,1 

— 0'k,2 

—<7fc,3 . . 

Cfc 


(A 3 ) 


< 7 i 2 appears in the four upper-left entries, including inside ci and C2, but nowhere else in 
the rest of the matrix. Using this fact and properties of determinants, we can rewrite the 
expression so that <712 will only appear in one place. Specifically we can replace the second 
row with the sum of itself and the first one and then repeat this procedure with the same 
columns 



Cl 

— <712 

— <7i3 


■<7l,fc 



c'l 

C'2 

— (<7i 3 + <723) ••• ~(<7l, 

k + <72,fc) 


det S' = 

—<731 

— <732 

C3 


■<73,fc 



— CTfc^l 

— 0'k,2 

~(^k,3 


Cfc 



Cl 


c'l 

— <7i3 

— 

<7l,fc 


c'l 

c'l 

T C2 — 

(<7i3 + <723) . 

• ~(<7i,A: + <72 


—O'si 

~(C31 + <732) 

C3 

— 



— <7fc,l 


1 + <7^,2) 

— <7fc,3 




(A 4 ) 
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Here = C\ — and C2 = C2 — cr^ do not contain cti 2. Hence the only entry that depends 
on (Ti 2 is Cl, and it does so linearly. We can expand the latter determinant as 


det S' = cri2 


c '2 -(0-13 +Cr23) 

■(0‘31+0'32) C 3 


-(^k,l + <^k,2) —^k,3 

+ (terms not containing (T12). 


— + cr2,k) 

—0'3,k 


(A 5 ) 


Dehne CT'2 := + crj^2 for i = 3 ,k. Then the determinant multiplying cri2 in Eq. (A 5 ) 


will be of the same form as the original determinant of S'. The size of this determinant is 
{k — 1 ) X [k — 1 ) (which corresponds to n = fc), so by the induction hypothesis it must be 
positive. Since the choice of ai2 was arbitrary we have proved that the coefficient in front 
of each edge conductance in S' is positive. Therefore S', as a polynomial in a’s, has only 
positive terms, which proves the lemma. □ 

Moreover, as S" has a form very similar to that of S', the proof above would work for its 
determinant as well. 


2. Equivalent conductance as a function of edge conductances 

From the consideration above it follows that both det S' and det S" are linear functions 
of edge conductances. Thus the equivalent conductance can be written as 

Acr^: + B 


<^eq — 


(A6) 


CcT* + D 

for an arbitrary edge conductance a*. Here H, B, C, and D depend on all the edge conduc¬ 
tances other than cr*. 

There is an important property of the equivalent conductance that is also worth men¬ 


tioning. In Sec. HC we saw that a special arrangement of some edge conductances ( 1 , 2 , 4 
and 5 ) resulted in zero current through the middle edge ( 3 ). In that case, it is easy to see 
that (Jeq does not depend on <73. Once such a symmetry is recognized, one can do two things 
without affecting the equivalent conductance: 

• Throw the middle edge away, i.e., set CT3 = 0 . This can be done, since there is no 
current through this edge. 
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• Short-circuit the top and bottom nodes, i.e., put as —)■ cx). This can be done, since 
the nodal potentials V2 and V3 are equal. 


In both cases, the resulting circuit can be easily simplihed and the equivalent conductance 
can be computed according to Eqs. and Q. Importantly, these two resulting circuits 
are different, but have the same equivalent conductance. This may not be as obvious for 
a more complicated circuit. Suppose there is special edge (with conductance a*) such that 
performing the two operations above yields the same equivalent conductance. We can prove 


then that (Jeq does not depend on a*, as follows. Requiring in (A6) that creq(O) = creq(oo) 
yields A/C = B/D, which implies that a* drops out from the equivalent conductance. 


Appendix B: Simplifiable circuits 

In this section we consider circuits that contain elements in series or in parallel, as well 
as circuits that can be reduced for symmetry reasons. Specihcally, we are interested in the 
form of the conductance matrix aij and the corresponding S-matrix. 


1. Connection in series 

If there is a pair of edges in series, the node {k) shared by these edges would be two-valent 
and the corresponding row/column in cr^ would have only two non-zero entries, say, aki and 
akm- This also implies that the nodes I and m are not connected directly, i.e., aim = 0 . From 
physics we know that the two edges can be replaced by one with the equivalent conductance 
given by ([^. Therefore, we can reduce the size of the conductance matrix by crossing out 
the row and column and by setting aim = <^mi '■= if^ki + 

It is also insightful to investigate this statement mathematically. Without loss of gener¬ 
ality, we can set k = 2 , I = 3 , and m = 4 . Then the (top-left part of the) S'-matrix for such 
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a circuit would look like 


^ Cl 0 —Uis —(Ti4 —(Ji5 .. 

0 C2 ~U23 —(724 0 ... 

— CTsi —(732 C3 0 —<735 . . . 

— C741 —(742 0 C4 —(745 • • • 

— (751 0 —(752 —(753 C5 

V : ! ! : ! 


(Bl) 


As explained above, (734 = (743 = 0. Since node 2 is only connected to 3 and 4 , the non-zero 
entries of the second row and column are (723 = (732, <J24 = (742 and C2 = (723 -|- (724. On the 
other hand, the matrix describing the reduced circuit is 


Cl 

— (7i3 

—(714 

— (7i5 • 

• • 1 


— C'31 

C3 

5'34 

— 0-35 . 



— (741 

^■43 

C4 

— (745 • 


(B 2 ) 

“C'51 

— CTss 

~cr54 

C5 . 




V ; : : : -.y 


Here (734 = ((733^ -|- (734^) ^ = c’'23(724/((723 -|- (724), whereas C3 and C4 include (734 = (743. We 
can similarly introduce the sub-matrices S" and S". We shall now prove the following 


Lemma 4 


det S' det S' 

“ det S" “ det S" ~ 


(B 3 ) 


Proof: In what follows, for the sake of compactness we omit the dots in the determinants. 


We start by expanding det S' from (Bl) in the elements of the second row 


Cl 

— (7i3 

— (7i4 

—(7i5 

— (731 

C3 

0 

—^35 

— (741 

0 

C4 

~Cr45 

“^■51 

— 0-53 

~cr54 

C5 



Cl 

0 

— (7i4 

— (7i5 


Cl 

— 7‘i3 

0 

—(7i5 


—cr3i 

— <732 

0 

— 0-35 

+ (724 

— C'31 

C3 

— <732 

—0-35 

+ (723 

— (741 

“^■42 

C4 

“^■45 

— (741 

0 

C4 

“<742 



— (751 

0 

— (754 

C5 


— (751 

— ^53 

0 

C5 


(B 4 ) 
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Note that we switched columns 2 and 3 in the last determinant. Using that C2 = CT23 + CT34, 
we can collect similar terms as 


(B6) 


It is now easy to see that the sum of columns 2 and 3 is the same for both determinants. 
Thus we can replace column 2 in each term with this sum, which would allow us to combine 
the two determinants into one 



Cl 

— <7 i 3 

—<7i4 

— <7 i 5 


Cl 

— <7 i 3 

— <7i4 

— <7 i 5 


—<731 

C3 — <732 

0 

— (735 


— <731 

C3 

— (732 

— (735 

det Y! = CT23 

—<741 

— <742 

C4 

— <745 

+ <724 

— <741 

0 

C4 — (J42 

— <745 


“<751 

— <753 

“<754 

C5 


~<751 

— (753 

— (754 

C5 


Cl 

-((713 + (714) 

— (<723 + <724)<7 i 4 

— <7 i 5 

— <731 

C3 — <723 

<723(724 

— <735 

— <741 

1 

to 

(<723 + <724) c 4 — < 7^4 

— <745 

— <751 

— (<753 + <754) 

— (<723 + <724)<754 

C5 


(B6) 


We can now factor ((T23 + cr24) out of the third column and use (CT23CT24)/(CT23 + CT24) = CT34 (in 
the second entry) together with o\J ((T23 + CT24) = CT24 — 5'34 (in the third entry) to obtain 


(<723 + <724) 


Cl 

-((7 i 3 + (714) 

— <7 i 4 

— <7 i 5 

— <731 

C3 — <723 

£■34 

— (735 

— <741 

C4 — CT42 

C4 — CT24 + £'34 

— <745 

— <751 

— (<753 + <754) 

— <754 

C5 


(B 7 ) 


where the third diagonal entry is precisely C4 — (T24 + a‘34 = £4. Similarly C3 — 02^ + cr34 = £3. 
With that in mind, subtracting the third column from the second one yields 


det Y! — (<723 + <724) 


Cl 

—(7 i 3 

— <7 i 4 

— <7 i 5 

— (731 

£3 

— £■34 

— (735 

— <741 

—£■34 

£4 

— <745 

— <751 

—(753 

— <754 

C5 


= (<723 + <724) det Yl. 


(B8) 


Repeating the same consideration without the hrst row and column we would get 


det Y!' — (<723 + <724) 


C3 

— <734 

—<735 

— £■34 

£4 

—<745 

— (753 

— <754 

C5 


= (<723 + <724) det Y!\ 


(B 9 ) 


from which (B 3 ) follows. 


□ 
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2. Connection in parallel 


If two nodes are connected by multiple edges, those edges can be replaced by one with 
the total conductance from ([^. This is the reason we did not consider multiple edges from 
the start, without any loss of generality. 

3. Short-circuiting two nodes 

The two above simplihcations do not affect the equivalent conductance of the whole 
circuit. A more interesting scenario (that does, in general) is short-circuiting two nodes, 
say i and j, by connecting them with an ideal wire. Mathematically this corresponds to 
letting (Tjj —)■ oo, whereas the physical implication is the equality of the corresponding nodal 
potentials V) = Vj. Such nodes can be simply merged together, becoming one. After the 
merging, some edges that were not parallel may become parallel. Thus we can use the idea 
of the previous paragraph. 

In essence, this short-circuiting eliminates one unknown potential, hence reduces the size 
of the conductance matrix. Instead of the two rows/columns {i and j) we have only one. 
The entries of this new row/column are given by 

= CTifc + o-jfc, (BIO) 


for any k ^ i,j. 

So far we have not discussed the diagonal elements of aij. One reason being is that they 
do not affect the E-matrix, which has been most relevant in this paper. In fact, we can 
generalize the procedure of going from to E as follows 


Eij = Diag(ci,C2,...,Cn) 


a. 




(Bll) 


with the same c’s as before: Ci= Yh ^ij- Clearly, no matter what an are, they drop out from 

i=i 

E, according to (Bll). 

At the same time, what would be the meaning of anl How can a node be connected 
to itself? Now, in the spirit of the short-circuiting scenario above, one can think of each 
node as “self short-circuited”. In other words, an = oo. In fact, this observation will 


make the short-circuiting recipe (BIO) valid for diagonal elements as well. Notice that all of 
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the manipulations with cxjj discussed in this appendix work in the same exact way for the 
S-matrix. In other words, these manipulations “commute” with (Bll). 


To conclude, we would like to comment of the symmetry issue. Some nodal potentials 
may turn out to be equal on symmetry grounds, even if the corresponding nodes are not 
short-circuited. For instance, as we pointed out in Section |V| the Wheatstone bridge circuit 


has V2 = V3, if ( 18 ) is satished. Another example is the classic resistor cube problem, where 
there are triples of (not connected) nodes having the same potential. Importantly, in these 
symmetric situations, the short-circuiting does not affect the equivalent conductance. As 
such symmetries are not as manifest in more complicated circuits, it would be interesting 
to come up with a way of detecting the equipotential nodes by looking at the form of the S 


matrix ( 33 ). 
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